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Abstract 

This article deals with the problem of optimal allocation of capital to corporate 
bonds in fixed income portfolios when there is the possibility of correlated defaults. 
Using a multivariate normal Copula function for the joint default probabilities we show 
that retaining the first few moments of the portfolio default loss distribution gives an 
extremely good approximation to the full solution of the asset allocation problem. We 
provide detailed results on the convergence of the moment expansion and explore how 
the optimal portfolio allocation depends on recovery fractions, level of diversification 
and investment time horizon. Numerous numerical illustrations exhibit the results for 
simple portfolios and utility functions. 



^To appear in, Correlated Default Analysis for Collateralized Debt Obligations, edited by S. Das, G. 
Fong and N. Kapadia. 



1 Introduction 

Investors routinely look to the corporate bond market in particular, and spread markets 
in general, to enhance the performance of their portfolios. However, for every source of 
excess return over the risk-free rate, there is a source of excess risk. When sources of 
risk are correlated, the allocation decision to the risky sectors, as well as allocation to 
particular securities in that sector, can be substantially different from the uncorrelated 
case. Since the joint probability distribution of returns of a set of defaultable bonds varies 
with the joint probabilities of default, recovery fraction for each bond, and the number of 
defaultable bonds, a direct approach to the allocation problem that incorporates all these 
factors completely can only be attempted numerically within the context of a default model. 
This approach can have the short-coming of hiding the intuition behind the asset allocation 
process in practice, which leans very heavily on the quantification of the first few moments, 
such as the mean, variance and skewness. In this paper, we will take the practical approach 
of characterizing the portfolio default loss distribution in terms of its moment expansion. 
Focusing on allocation to corporate bonds (the analysis in this paper can be generalized 
to any risky sector which has securities with discrete payoffs), we will answer the following 
questions: 

• In the presence of correlated defaults, how well does retaining the first few moments 
of the portfolio default loss distribution do, for the portfolio allocation problem, as 
compared to the more intensive full numerical solution? 

• For different choices of correlations, probabilities of default, number of bonds in the 
portfolio, and investment time horizon, how does the optimal allocation to the risky 
bonds vary? 

In this paper we consider portfolios consisting of risk free assets at a return y and 
corporate zero coupon bonds with return Cj for firm i, and study how correlations between 
defaults affect the optimal allocation to corporate assets in the portfolio over some time 
horizon T. We assume that if company i defaults at time t a fraction Ri (the recovery 
fraction) of the value of its bonds is recovered and reinvested at the risk free rate. We will 
quantify the impact of risk associated with losses from corporate defaults ^ on portfolio 
allocation. To compensate the investor for default risk, the corporate bond's return q is 
greater than that of the risk-free assets in the portfolio. The excess return, or spread of 
investment in bonds of firm i can be decomposed into two parts. The first part of the spread, 
which we call Aj, arises as the actuarially fair value of assuming the risk of default. The 
second part, which we call /Xj, is the excess risk premium that compensates the corporate 
bond investor above and beyond the probabilistically fair value. In practice, can arise due 
to a number of features not directly related to defaults. For example, traders will partition 
/ij into two pieces, one for liquidity and one for event risk, ^Uj = /j + e,. Most low grade 
bonds may have a full percentage point arising simply from liquidity premia, and /ij can 
fluctuate to large values in periods of credit stress. Liquidity li is systematic and one would 
expect it to be roughly equal for a similar class of bonds. Event risk contributes to due 

^In this paper, we use annualized units so quoted default probabilities, hazard rates and corporate rates 
are annualized ones. 
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to the firm's specific vulnerability to factors that affect it (e.g. negative press). In periods 
of stress the default probabilities, li and all increase simultaneously, and the recovery rate 
expectations Ri fall, leading to a spike in the overall spread. Since these variables can be 
highly volatile, the reason behind the portfolio approach to managing credit is to minimize 
the impact of non-systematic event risk in the portfolio.^ The excess risk premium itself is 
not very stable over time. Empirical research shows that the excess risk premium might vary 
from tens of basis points to hundreds of basis points. For instance, in the BB asset class, if 
we assume a recovery rate of 50% and default probability of 2%, the actuarially fair value 
of the spread is 100 basis points. However, it is not uncommon to find actual spreads of the 
BB class to be 300 bp over treasuries [Altman (1989)]. The excess 200 bp of risk premium 
can be decomposed in any combination of liquidity premium and event risk premium, and 
is best left to the judgment of the market. When the liquidity premium component is small 
compared to event risk premium, we would expect that portfolio diversification and the 
methods of this paper are valuable. 

One other factor needs to be kept in mind when comparing historical spreads to current 
levels. In the late eighties and early nineties, the spread was routinely quoted in terms of 
a treasury benchmark curve. However, the market itself has developed to a point where 
spreads are quoted over both the libor swap rate and the treasury rate, and the swap rate 
has gradually substituted the treasury rate as the risk-free benchmark curve. This has two 
impacts. Firstly, since the swap spread (swap rate minus treasury rate for a given maturity) 
in US is significant (of the order of 50 bp as of this writing), the excess spread needs to 
be computed as a difference to the swap yield curve. Secondly, the swap-spread itself has 
been very volatile during the last few years, which leads to an added source of non default 
related risk in the spread of corporate bonds when computed against the treasury curve. 
Thus, the 200 bp of residual spread is effectively 150 bp over the swap rate in the BB 
example, of which, for lack of better knowledge, equal amounts may be assumed to arise 
from liquidity and event risk premium over the long term. The allocation decision to risky 
bonds strongly depends on the level of risk-aversion in the investor's utility function, and 
the required spread for a given allocation will go up nonlinearly as risk aversion increases. 

The value of the optimal fraction of the portfolio in corporate bonds, here called aopt, 
cannot be determined without knowing the excess risk premium part of the corporate re- 
turns. They provide the incentive for a risk averse investor to choose corporate securities 
over risk free assets. In this paper we explore, using utility functions with constant relative 
risk aversion, the convergence of the moment expansion for Oopt- Our work indicates that, 
(for Oopt less than unity) Oopt is usually determined by the mean, variance and skewness 
of the portfolio default loss probability distribution. The sensitivity to higher moments 
increases as Uopt does. Some measures of default risk, for example a VAR analysis^, may be 
more sensitive to the tail of the default loss distribution. We also examine how the optimal 
portfolio allocation scales with the number of firms, time horizon and recovery fractions. 

Historical evidence suggests that on average default correlations increase with the time 
horizon^. For example, Lucas (1995) estimates that over one year, two year and five year 

^The authors would like to thank David Hinman of PIMCO for enlightening discussions on this topic. 
^See for example, Jorion (2001). 

^ Zhou (2001) derives an analytic formula for default correlations in a first passage-time default model 
and finds a similar increase with time horizon. 
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time horizons default correlations between Ba rated firms are 2%, 6% and 15% respectively. 
However, the errors in extracting default correlations from historical data are likely to be 
large since defaults are rare. Also these historical analysis neglect firm specific effects 
that may be very important for portfolios weighted towards a particular economic sector. 
Furthermore, in periods of market stress default probabilities and their correlations increase 
[Das, Freed, Geng and Kapadia (2001)] dramatically. 

There are other sources of risk associated with corporate securities. For example, the 
market's perception of firm i's probability of default could increase over the time horizon T 
resulting in a reduction in the value of its bonds. For a recent discussion on portfolio risk 
due to downgrade fluctuations see Dynkin, Hyman and Konstantinovsky (2002). Here we 
do not address the issue of risk associated with fluctuations in the credit spread but rather 
focus on the risk associated with losses from actual defaults. 

In the next section a simple model for default is introduced. The model assumes a 
multivariate normal Copula function for the joint default probabilities. In section 3 we 
set up the portfolio problem. Moments of the fractional corporate default loss probability 
distribution are expressed in terms of joint default probabilities and it is shown how these 
can be used to determine Oopt- In section 4 the impact of correlations on the portfolio 
allocation problem is studied using sample portfolios where all the firms have the same 
probabilities of default and the correlations between firms are all the same. The recovery 
fractions are assumed to be zero for the portfolios in section 4. The impact of non-zero 
recovery fractions on the convergence of the moment expansion is studied in section 5. 
Concluding remarks are given in section 6. 

This work is based on Wise and Bhansali (2002). It extends the results presented in 
that paper to arbitrary time horizons (allowing default to occur at any time) and makes 
more realistic assumptions for the consequences of default. 



2 A Model For Default 

It is convenient for discussions of default risk to introduce the random variables ni(tj). 
ni{ti) takes the value 1 if firm i defaults in the time horizon ti and zero otherwise. The 
joint default probabilities are expectations of products of these random variables, 

Ph..A,^{th,---:ti^)=E[hi^{ti^)---hi^{tiJ], when / is • • • / (1) 

Pi{ti) is the probability that firm i defaults in the time period ti and Pii,...,i^(ti^, ■ ■ ■ ,ti^) 
is the joint probability that the m-firms ii, . . Am default in the times periods ti-^, . . . ,ti^. 

We assume that the joint default probabilities are given by a multivariate normal Copula 
function [See for example, Lee (2000)]. Explicitly, 



1 r-xiih) r-Xn{tn) 

Pi...„(ti, . . . ,tn) = - — -^-7==/ dxi--- dxnexp 



(27r)2 VdeU-'-oo 



00 



1 
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2 



(2) 

where the sum goes over i,j = 1, . . . ,n, and is the inverse of the n x n correlation 

matrix ^ij. For n not too large the integrals in equation (§) can be done numerically or, since 
defaults are rare, analytic results can be obtained using the leading terms in an asymptotic 
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expansion of the integrals. The choice of a multivariate normal Copula function is common 
but somewhat arbitrary. In principle the Copula function should be chosen based on a 
comparison with data. For recent work along these lines see Das and Ceng (2002). 

In equation ^ the nxn correlation matrix is usually taken to be the asset correlation 
matrix. This has the advantage of allowing a connection to stock prices [Merton (1974)]. 
However this assumption is not necessary. For example the assets, aj could be functions 
of normal random variables, Oj = gi{zi). Suppose default occurs if the assets Oj cross the 
thresholds Tj. Then the condition for default on the normal production factor variables Zi 
is that they cross the thresholds g\ ^\Ti). In such a model it is natural to interpret 
as the correlation matrix for the normal production factor variables Zj. If the functions gi 
are linear then the assets are also normal and their correlation matrix is also given by ^ij. 
However if the functions gi are not linear the joint probability distribution for the assets is 
not multivariate normal and can have fat tails. Unless the gi are specified it is not possible 
to connect stock prices to the correlation matrix and the default thresholds Xi- However, 
even when the functions gi are not known the default thresholds or "equivalent distances 
to default" XiiU) and the correlation matrix are determined by the default probabilities 
Pi{ti) and the default correlations dij[t) and so have a direct connection to measures that 
investors use in quantifying security risk. In our work we assume that the correlation matrix 
is time independent. 

Equation (pf) in the case n = 1 gives. 



Pi{ti) 



1 



-XiiU) 



dxi exp 



1 



(3) 



Hence for an explicit choice for the time dependence of the default probabilities the default 
thresholds are known. A simple (and frequently used) choice for the default probabilities is 



Pi{ti) = 1 - exp(- 



(4) 



where the hazard rates hi are independent of time. 

Since nj(tj)^ = ni[ti) it follows that the correlation of defaults between two different 
firms (which we choose to label 1 and 2) is, 



E[hi{t)h2{t)] - E[ni{t)]E[n2{t)] 



Pl2it,t) - Pi{t)P2{t) 

VPimi-Pi{tmmi-P2{t)) 



(5) 



The default model we are adopting is not as well motivated as a first passage-time default 
model [Black and Cox (1976), Longstaff and Schwartz (1995), Leland and Toft (1996), etc.] 
where the assets undergo a random walk and default is associated with the first time that 
the assets fall below the liabilities. However, it is very convenient to work with. 

We will evaluate the integrals in equation (|2|) by numerical integration. For this we 
need the inverse and determinant of the correlation matrix. It is very important that the 
correlation matrix is positive semi-definite. If it has negative eigenvalues the integrals 
in equations @ are not well defined. Typically a correlation matrix that is forecast using 
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qualitative methods will not be mathematically consistent and have some negative eigen- 
values. A practical method for constructing the consistent correlation matrix that is closest 
to a forecasted one is given in Rebonato and Jackel (2000). For the portfolios discussed 
in sections 4 and 5 the n x n correlation matrix ^ij is taken to have all of its off diagonal 
elements the same, i.e., = ^ for i ^ j. Such a correlation matrix has one eigenvalue 
equal to 1 + (n — 1),^ and the others equal to 1 — ^. Consequently its determinant is 

det[^,,-] = (l-er"^[l + (n-l)C]. (6) 
Its inverse has diagonal elements, 

^ l + {n-2)^ 

l + (n-2)C-(n-l)e2' ^ ' 

and off diagonal elements {i ^ j), 

^(71) = i (8) 

1 + (n - 2)e - (n - 1)^2 ^ ' 

In the next section we consider the problem of portfolio allocation for portfolios con- 
sisting of corporate bonds subject to default risk and risk free assets. The implications of 
default risk are addressed using the model discussed in this section. Other sources of risk, 
for example, systematic risk associated with the liquidity part of the excess risk premium, 
are neglected. 



3 The Portfolio Problem 

Assume a zero coupon bond from company i grows in value (if it doesn't default) at the 
(short) corporate rate Cj(t) and that if the company defaults a fraction Ri of the value of 
that bond at the time of default is reinvested at the risk free (short) rate y{t). Then, the 
random variable for the value of this bond at some time T in the future is, 

|2 = exp^"<irc.M)(l-%(r)) 

+ Ri£ds^^exp(^l^drCi(T)dT + £dTy{T)\, (9) 

where Vi{0) is the initial value of the zero coupon bond. In equation (P) we have continuously 
compounded the returns and we have assumed that T is less than or equal to the maturity 
date of the bond^. Note that the random variable dhi(t)/dt is equal to the Dirac delta 
function, 5(t — ti), if company i defaults at the time ti and is zero if it doesn't default. The 
first term in equation gives the value if the company does not default before time T and 

^Assuming that after maturity the value of a zero coupon corporate bond is reinvested at the risk free 
rate it is straightforward to generaUze the analysis of this paper to portfolios containing zero coupon bonds 
with different maturities, some of which are less than the investment horizon. Similarly default risk for 
portfolios containing coupon paying bonds can be studied since each coupon payment can be viewed as a 
zero coupon bond. 



5 



the second term (proportional to Ri) gives the value if the company defaults before time T. 
The corresponding formula for a risk free asset is, 



VrfjT) 
VrfiO) 



exp 



(10) 



Note that we are not treating the (short) risk free rate y(r) as a random variable, although 
it is certainly possible to generalize this formalism to do that. The (short) corporate rate 
is decomposed as 

Ci{t)=y{t) + Xi{t)+fii{t), (11) 

where /ij is the excess risk premium and Aj is the part of corporate short rate that compen- 
sates the investor in an actuarially fair way for the fact that the value of the investment 
can be reduced through corporate default. In other words, 



E[Vi{T)/Vim\,^=0 = Vrf{T)/Vrf{0). 

Taking the expected value equation (|l^) implies that 

dPiis) 



1 = (1 - Pi{T)) exp 1^^^ dTA,(r)^ + Ri ds- 



ds 



exp 



drXdr] 



Differentiation with respect to T gives 



dP,{T) l-Ri 



(12) 



(13) 



(14) 



dT l-Pi{T)- 

for the time dependence of the probabilities of default equation (|14| 

R,.). 



Using equation 
implies the familiar relation 

Xi{T) = h,{l-R,). (15) 

The above results imply that the total wealth, after time T, in a portfolio consisting of 
risk free assets and zero coupon corporate bonds that are subject to default risk is (assuming 
the bonds mature at dates greater than or equal to T) 



W{T) 



Wq exp 



dryir) x 



(1 - a) + a (^fisxp (^J^ 



dr(A,(r)+^,(T)) -1{T) 



(16) 



where, Wq is the initial wealth, a is the fraction of corporate assets in the portfolio, 1{T) is 
the random variable for the fractional default loss over the time period T, Ri is the recovery 
fraction and fi denotes the initial fraction of corporate assets in the portfolio that are in 
firm i (J2i fi — !)• equation (16) the sum over i goes over all N firms in the portfolio 
and the default loss random variable is given by, 



f(T) 



i 

R. 



exp 







^ dhi{s) 
ds — ; — exp 
ds 







dr(A,(r) + ^,(T))j n,(T) 

dr((Ai(T) +^i(r)) 



(17) 
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Integrating by parts, 

i{T) = E/^ 



T 



exp I / (ir(A,(r) + /ii(r))j ni{T){l - R,) 
+ Ri dshi{s){Xi{s) + iJi{s))ex.p (^J^ dT{{Xi{T) + fii{T)) 
The expected value of 1{T) is, 



(18) 



mm] = Y.f^ 



exp / dr(A,(r)+;u,(T)) Pi(r) 



T 

Rj I ds 



dPjjs ) 
ds 



exp / dT{{Xi{T) + ^i(T)) 



(19) 



Fluctuations of the fractional loss l(T) about its average value are described by the variable, 

61{T) = i{T) - E[l{T)]. (20) 

The mean of 61 is zero and the probability distribution for 61 determines the default risk of 
the portfolio associated with fluctuations of the random variables hi{t). The moments of 
this probability distribution are, 



yi^'^T) = E[{6l{T)y 



(21) 



Using equation (|l|) and the property ni{ti)ni{t2) = nj(min[ti, t2])) the moments v^"^\T) 
can be expressed in terms of the joint default probabilities. 
The random variable for the portfolio wealth at time T is. 



T 



W{T) = Wo exp / dryir) [1 + ax(r) - a6l{T)], 



where. 



x(T) = -l + E/iexp dT{X,{T) + f,,{T))^ - E[liT)]. 
Using equations (|l3|) and (|l9|) we find that, 

x(T) = Y.fi 



(22) 



(23) 



exp 



dTfli{T)^ - 1 



+ Ri ds 



T 



dFiis) 
ds 



exp 



dr^Xiir) + /ii(r)) 1 - exp / drfnir) 



.(24) 



Note that x(T) vanishes if the excess risk premiums /Xj do. 

To find out what value of a is optimal for some time horizon T a utility function 
is introduced which characterizes the investor's level of risk aversion. Here we use utility 
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functions of the type U^{W) = /j which have constant relative risk aversio n[], 1-7. The 
optimal fraction of corporates, Oopt, maximizes the expected utility of wealth E[U^{W{T))]. 
Expanding the utility of wealth in a power series in 61 and taking the expected value gives 



E[U,{W{T))] 



(W^/j) exp 7 / dryir) (1 + ax(r))' 



V- r(m - 7) 
^„ r(-7)r(m + 1) V 1 + ax(r) 



m=2 



(25) 



where T is the Euler Gamma function. The approximate optimal value of a obtained from 
truncating the sum in equation (p5|) at the m'th moment is denoted by Um- The focus of 
this paper is on portfolios that are not leveraged and have Oopt less than unity. We will 
later see that typically, for such portfolios, the Om converge very quickly to Oopt so that for 
practical purposes only a few of the moments v^^^ need be calculated. Note that the value 
of Gopt is independent of the risk free rate y{t). 

The expressions for 1{T) and x(r) are complicated but they do simplify for short times 
T or if all the recovery fractions are zero. For example, if T is small 



z(r)~^/i(i-i?,)ni(r). 



(26) 



and 



x(r)=.^/, r dsi,,{s). 

, Jo 



On the other hand if all the recovery fractions vanish (i.e. Ri = 0) then, 



dr(A,(T)+/ii(r)) hi{T), 



and 



x(T) = J2h 



exp 



1 



(27) 



(28) 



(29) 



If x(T) is zero then aopt is also zero. It is x(r) that contains the dependence on the excess 
risk premiums which provide the incentive for a risk adverse investor to choose corporate 
bonds over risk free assets. An approximate formula for aopt can be derived by expanding 
it in x(T), 



a 



opt 



y £!i^px(^)^ 

^ n! 



n=l 



(30) 



The coefficients Si{T) can be expressed in terms of the moments, v^'"^\T). Explicitly, for 
the first two coefficients. 



(l-7)t;(2)(r)^ 



(31) 



^See, for example, IngersoU (1987). 
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and 

(2-7y3)(r) 
^^^^) = -(i-7)2(.(2)(r))3- 

For 7 < 1 including the third moment reduces the optimal fraction of corporates when 
v^^\T) is positive. 



4 Sample Portfolios with Zero Recovery Fractions 

Here we consider very simple sample portfolios where the correlations, default thresholds and 
excess risk premiums are the same for all N firms in the portfolio i.e., S^ij = ^, Xi(^) = x{t) 
and = Then all the probabilities of default are the same, Pi{t) = P{t), and 

the joint default probabilities are also independent of which firms are being considered, 
• • • tm) = -Pi2...m(ii) • • • ^m)- The time dependence of the default probabilities is 
taken to be given by equation (^. Given our assumptions the hazard rates are then also 
independent of firm, hi = h. We also take the portfolios to contain equal assets in the 
firms so that, fi = 1/N, for all i, and assume all the recovery fractions are zero. Since the 
portfolio allocation problem does not depend on the risk free rate y or the initial wealth we 
set y = and Wq = 1- 

With these assumptions the expressions for 1{T) and x(T) become, 

1{T) = ^ exp{hT + ^T) J2 f^iiT), (33) 

i 

and 

x(r) = exp(;ur) - 1. (34) 

For random defaults the probability of a fractional loss of I = exp(hT + fj,T)n/N in 
the time horizon T is (1 — P(T))'^~'^P(T)"'Nl/{N — n)ln\ and so the m'th moment of the 
default loss distribution is 



(-)(r) = exp(m/ir + m/.r)^(--p(r)j (i _ P(r))^-P(r)"^^^— ^. (35) 



n=0 



The expected utility of wealth is 

TV 



E[u,iwiT))] = i^(i-p(r))^-"p(r)" 

1 [N-ny.n] 
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1 + ax(r) -ai^—- P{T)j expihT + ^T) j . (36) 

Having the explicit expression for the utility of wealth in equation (|3^) lets us compare 
results of the moment expansion for the optimal fraction of corporates a^n with the all 
orders result, Oopt- These are presented in Table I in the case h = 0.02, fi = lOObp and 

*We also assume that fi is independent of time. 
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7 = —4. Results for different values of the number of firms and different time horizons 
T are shown in Table I. We also give in columns three and four of Table I the volatility 

vol = ^Jv(^){T), (37) 

and skewness 

skew = ^ (38) 

of the portfolio default loss probability distribution. 

Increasing gives a larger value for the optimal fraction of corporates because diversi- 
fication reduces risk. This occurs very rapidly with N. For N = 1, ^ = lOObp, T = lyr and 
7 = —4 the optimal fraction of corporates is only 7.7%. By = 10 increasing the number 
of firms has reduced the portfolio default risk so much that the lOObp excess risk premium 
causes a portfolio that is 76% corporates to be preferred. 

For all the entries in Table I the moment expansion converges very rapidly, although for 
low N it is the small value of aopt that is driving the convergence. Since in equation ( p5|) the 
term proportional to v^"^^ has a factor of a"* accompanying it we expect good convergence 
of the moment expansion at small a. 

The focus of this paper is on portfolios that are not leveraged and have Oopt < 1- But a 
value Oopt > 1 is not forbidden when finding the maximum of the expected utility of wealth. 
This occurs at lowest order in the moment expansion in the last row of Table I. 



Table I: Optimal Fraction of Corporates for Sample Portfolios with Random (^ = 0) 
Defaults. Other parameters used are: h = 0.02, 7 = —4 and /i = lOObp. 



T 


N 


vol 


skew 


02 


/ lyr 


1 


0.14 


6.9 


0.098 


lyr 


5 


0.064 


3.1 


0.50 


lyr 


10 


0.045 


2.2 


1.0 


5yr 


1 


0.34 


2.8 


0.090 


Syr 


5 


0.15 


1.2 


0.49 


\5yr 


10 


0.11 


0.87 


1.1 



"3 


04 


05 


Oiopt 


0.079 


0.077 


0.077 


0.077\ 


0.40 


0.39 


0.38 


0.38 


0.81 


0.78 


0.77 


0.76 


0.074 


0.072 


0.072 


0.072 


0.39 


0.37 


0.36 


0.36 


0.85 


0.75 


0.73 


0.71 / 



Next we consider the more typical case where there are correlations, ^ 7^ 0. Table Ila 
gives values of 02 to 05 for T = lyr, ^ = lOObp, h = 0.02 and 7 = —4. Portfolios with 
N = 10, 50 and 100 are considered and values for ^ between 0.5 and 0.25 are used. 

Again the convergence of the moment expansion is quite good. Although just including 
the variance {i.e. 02) can be off by almost a factor of two 03 is usually a reasonable 
approximation to the true optimal fraction of corporates. The convergence is worse the 
larger the value of aopt and for values of aopt less than 40% we find that as is within about 
10% of aopt- 

In the last row of Table Ila the value of 05 is equal to unity. However, we know that the 
true value of Oopt must be less than unity. For a = 1 there is some finite (but tiny) chance 
of the investor loosing all his wealth and for 7 < the utility of zero wealth is —00. 



10 



Table Ila: Moment Expansion for Optimal Portfolio with Correlated Defaults for T = lyr, 

/X = lOObp, h = 0.02 and 7 = -4. 



1 




IJ' 


? 


7 


1 

vol 


skew 


a2 


0:3 


a4 


"5 


/ lyr 


10 


lOObp 


0.50 


0.152 


0.070 


5.2 


0.42 


0.31 


0.29 


0.29 \ 


lyr 


10 


lOObp 


0.45 


0.126 


0.066 


4.9 


0.47 


0.35 


0.33 


0.32 


lyr 


10 


lOObp 


0.40 


0.102 


0.063 


4.5 


0.52 


0.39 


0.36 


0.36 


lyr 


10 


lOObp 


0.35 


0.082 


0.060 


4.2 


0.57 


0.43 


0.40 


0.40 


lyr 


10 


lOObp 


0.30 


0.064 


0.057 


3.8 


0.63 


0.48 


0.45 


0.44 


lyr 


10 


lOObp 


0.25 


0.049 


0.054 


3.5 


0.70 


0.53 


0.50 


0.49 


lyr 


50 


lOObp 


0.50 


0.152 


0.059 


5.5 


0.59 


0.42 


0.38 


0.37 


lyr 


50 


lOObp 


0.45 


0.126 


0.054 


5.2 


0.70 


0.49 


0.45 


0.43 


lyr 


50 


lOObp 


0.40 


0.102 


0.050 


4.8 


0.84 


0.58 


0.53 


0.51 


lyr 


50 


lOObp 


0.35 


0.082 


0.045 


4.4 


1.0 


0.70 


0.63 


0.61 


lyr 


50 


lOObp 


0.30 


0.064 


0.041 


4.0 


1.2 


0.85 


0.77 


0.74 


lyr 


50 


lOObp 


0.25 


0.049 


0.037 


3.6 


1.5 


1.1 


0.94 


0.90 


lyr 


100 


lOObp 


0.50 


0.152 


0.057 


5.6 


0.62 


0.44 


0.40 


0.39 


lyr 


100 


lOObp 


0.45 


0.126 


0.053 


5.3 


0.75 


0.52 


0.47 


0.45 


lyr 


100 


lOObp 


0.40 


0.102 


0.048 


5.0 


0.91 


0.62 


0.56 


0.54 


lyr 


100 


lOObp 


0.35 


0.082 


0.043 


4.6 


1.1 


0.76 


0.68 


0.65 


lyr 


100 


lOObp 


0.30 


0.064 


0.039 


4.2 


1.4 


0.94 


0.84 


0.80 


Vlyr 


100 


lOObp 


0.25 


0.049 


0.035 


3.8 


1.8 


1.2 


1.1 


1.0 J 



Correlations dramatically effect the dependence of the optimal portfolio allocation on 
the total number of firms. For the ^ = 0.5 entries in Table Ila the optimal fraction of 
corporates is, 0.29, 0.37, and 0.39 for N = 10, 50 and 100 respectively. For N = 10, 000 
we find that 05 = 0.40. Increasing the number of firms beyond 100 only results in a small 
increase in the optimal fraction of corporates. When defaults are random the moments of 
the portfolio default loss distribution go to zero as iV — 00. For example, with ^ = the 
variance of the default loss distribution is, 

v^'^^ = eM^hT + 2fiT) ^^^^ , (39) 
and the skewness of the default loss distribution is 

^;(3) 1 

VNP{l-P)^^~^^^' ^^^^ 

For random defaults as ^ 00 the distribution for 51 approaches the trivial one where 
SI = occurs with unit probability. However, for ^ > the moments v^"^^ go to non- 
zero values in the limit N ^ 00 and the default loss distribution remains non-trivial and 
non-normal. 

In Table lib the time horizon is changed to T = 5yr but the other parameters are left 
the same as in Table Ila. The value of aopt is always smaller than in Table Ila indicating 
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that the compounding of the excess risk premium does not completely compensate for 
the added risk associated with the greater probability of default. The convergence of the 
moment expansion and the dependence of the optimal fraction of corporates on the number 
of bonds is similar to that in Table Ila. 

Table lib: Moment Expansion for Optimal Portfolio with Correlated Defaults for T = 5yr, 

/i = lOObp, h = 0.02 and 7 = -4. 



T 


N 






dij 


vol 


skew 


02 


as 


a4 


as 


/ 5yr 


10 


lOObp 


0.50 


0.246 


0.193 


2.3 


0.29 


0.22 


0.21 


0.20 \ 


5yr 


10 


lOObp 


0.4S 


0.213 


0.184 


2.2 


0.32 


0.24 


0.23 


0.22 


5yr 


10 


lOObp 


0.40 


0.182 


0.175 


2.1 


0.36 


0.27 


0.25 


0.25 


5yr 


10 


lOObp 


0.3S 


0.153 


0.166 


2.0 


0.40 


0.30 


0.28 


0.27 


5yr 


10 


lOObp 


0.30 


0.127 


0.158 


1.9 


0.45 


0.34 


0.31 


0.31 


Syr 


10 


lOObp 


0.2s 


0.102 


0.149 


1.8 


0.51 


0.38 


0.35 


0.34 


5yr 


50 


lOObp 


0.50 


0.246 


0.174 


2.4 


0.36 


0.27 


0.25 


0.24 


Syr 


50 


lOObp 


0.4S 


0.213 


0.163 


2.3 


0.42 


0.31 


0.28 


0.28 


Syr 


50 


lOObp 


0.40 


0.182 


0.152 


2.2 


0.49 


0.35 


0.32 


0.31 


Syr 


50 


lOObp 


0.3S 


0.153 


0.141 


2.1 


0.58 


0.41 


0.38 


0.36 


Syr 


50 


lOObp 


0.30 


0.127 


0.129 


2.0 


0.71 


0.49 


0.44 


0.43 


Syr 


50 


lOObp 


0.2s 


0.102 


0.118 


1.8 


0.89 


0.60 


0.S4 


0.51 


Syr 


100 


lOObp 


0.50 


0.246 


0.172 


2.4 


0.38 


0.28 


0.26 


0.2s 


Syr 


100 


lOObp 


0.4S 


0.213 


0.160 


2.3 


0.44 


0.32 


0.29 


0.28 


Syr 


100 


lOObp 


0.40 


0.182 


0.149 


2.2 


0.52 


0.37 


0.33 


0.33 


Syr 


100 


lOObp 


0.3S 


0.153 


0.137 


2.1 


0.62 


0.43 


0.39 


0.38 


Syr 


100 


lOObp 


0.30 


0.127 


0.12s 


2.0 


0.76 


0.S2 


0.47 


0.4S 


\Syr 


100 


lOObp 


0.2s 


0.102 


0.113 


1.8 


0.98 


0.65 


0.S7 


0.54/ 



Table lib indicates that the volatility of the portfolio default loss probability distribution 
increases as the time horizon increases and the skewness decreases as the time horizon in- 
creases. The same behavior occurs for random defaults and this can be seem from equations 
(|39| ) and (|4^) with the small probability of default P increasing with time. 

Next we consider the effect of changing the utility function to one with a lower level 
of risk aversion. In Table III 7 = — 2 is used. Equation (^) suggests that an excess 
risk premium of 60bp would yield roughly the same portfolio allocation as the excess risk 
premium of lOObp used in Tables II. In Table III we consider portfolios with 100 firms and 
take the same parameters as in Table Ila apart from the lower level of risk aversion, 7 = —2 
and the lower excess risk premium, fi = 60bp. 

It is not surprising that the values of a2 in Table III are close to those in Table Ila. 
However, 05 is typically about 10% larger than in Table Ila. The convergence of the moment 
expansion is better in Table III than in Table Ila. 
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Table III: Moment Expansion for Optimal Portfolio with Correlated Defaults for T = lyr, 

N = 100, /X = 60bp, h = 0.02, and 7 = -2. 



T 


N 






dij 


vol 


skew 


0(2 


as 


a4 


05 




100 


60bp 


0.50 


0.152 


0.057 


5.6 


0.62 


0.47 


0.45 


0.44 \ 


lyr 


100 


60bp 


0.45 


0.126 


0.053 


5.3 


0.74 


0.56 


0.53 


0.52 


lyr 


100 


60bp 


0.40 


0.102 


0.048 


5.0 


0.90 


0.68 


0.63 


0.62 


lyr 


100 


60bp 


0.35 


0.082 


0.043 


4.6 


1.1 


0.82 


0.76 


0.74 


lyr 


100 


60bp 


0.30 


0.064 


0.039 


4.2 


1.4 


1.0 


0.94 


0.92 


Viyr 


100 


60bp 


0.25 


0.049 


0.035 


3.8 


1.7 


1.3 


1.2 


1.1 / 



The examples in this section have all used an annualized hazard rate of 2%, however, the 
convergence of the moment expansion is similar for significantly larger default probabilities. 
In Table IV the convergence of the moment expansion is examined for portfolios with 100 
corporate bonds and hazard rates of 2%, 4%, 6%, 8%, and 10%. The value of ^ is fixed at 
0.50 and the time horizon is lyr. The convergence of the moment expansion is good for all 
the values of h used in Table IV. 

For longer time horizons the increased default risk is not compensated for by the com- 
pounding of the excess risk premium. For example if we take the parameters of the last row 
of Table IV but increase the time horizon to T = 5yr then, a2 = 0.25, as = 0.21, = 0.18 
and a5 = 0.17. 

Table IV: Probability of Default and Convergence of the Moment Expansion for T = lyr, 

( = 0.50, N = 100 and 7 = -4. 



T 


h 






dij 


vol 


skew 


a2 


as 


a4 


"5 




0.02 


lOObp 


0.50 


0.152 


0.057 


5.6 


0.62 


0.44 


0.40 


0.39 \ 


lyr 


0.04 


200bp 


0.50 


0.189 


0.092 


4.0 


0.50 


0.36 


0.33 


0.32 


lyr 


0.06 


300bp 


0.50 


0.214 


0.121 


3.2 


0.44 


0.32 


0.29 


0.29 


lyr 


0.08 


500bp 


0.50 


0.231 


0.148 


2.7 


0.52 


0.35 


0.32 


0.31 


Viyr 


0.10 


600bp 


0.50 


0.246 


0.173 


2.4 


0.46 


0.32 


0.29 


0.28 / 



5 Sample Portfolios with Non-Zero Recovery Fractions 

In this section we consider portfolios like those in section 4 but allow the recovery fractions 
Ri = Rio be nonzero. Here we focus on short time horizons using the approximate formulas, 

iiT) = {l-R)hi{T), (41) 

and 

x(T) = ^T, (42) 

appropriate to that regime. Recall that if R is not zero the actuarially fair corporate spread 
is A = h{l — R). Table V shows the impact of nonzero recovery on agpt- The parameters 
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used are 7 = -4, T = lyr, N = W, R = 0.5, h = 0.04 and = 50 and lOObp. Doubling the 
hazard rate to 0.04 gives an actuarially fair corporate spread of 200bp, which is the same 
as in Tables II. Note that even though the actuarially fair spread in Table V is the same as 
in Table Ila nonzero recovery reduces the default risk and so the excess risk premium must 
be reduced by a factor of two to get values of Ogpt close to those in Table Ila. 

Table V: Recovery and the Moment Expansion for Optimal Portfolio with Correlated 



Default 




T = lyr, N = 




04 


R = 0.5 and 


7 


-4. 


T 


N 






dij 


vol 


skew 


02 


03 




"5 




10 


50bp 


0.50 


0.152 


0.050 


4.6 


0.40 


0.32 




0.31 \ 


lyr 


10 


50bp 


0.45 


0.126 


0.048 


4.5 


0.44 


0.36 


0.34 


0.34 


lyr 


10 


50bp 


0.40 


0.102 


0.046 


4.3 


0.49 


0.39 


0.38 


0.38 


lyr 


10 


50bp 


0.35 


0.082 


0.043 


4.2 


0.54 


0.44 


0.42 


0.42 


lyr 


10 


50bp 


0.30 


0.064 


0.041 


4.0 


0.60 


0.48 


0.47 


0.46 


lyr 


10 


50bp 


0.25 


0.049 


0.039 


3.8 


0.67 


0.54 


0.52 


0.51 


lyr 


10 


lOObp 


0.50 


0.152 


0.050 


4.6 


0.81 


0.57 


0.53 


0.52 


lyr 


10 


lOObp 


0.45 


0.126 


0.048 


4.5 


0.90 


0.63 


0.59 


0.57 


lyr 


10 


lOObp 


0.40 


0.102 


0.046 


4.3 


1.0 


0.70 


0.65 


0.63 


lyr 


10 


lOObp 


0.35 


0.082 


0.043 


4.2 


1.1 


0.78 


0.72 


0.70 


lyr 


10 


lOObp 


0.30 


0.064 


0.041 


4.0 


1.2 


0.86 


0.80 


0.78 


\lyr 


10 


lOObp 


0.25 


0.049 


0.039 


3.8 


1.4 


0.96 


0.88 


0.86 / 



For random defaults one can easily see why recovery reduces portfolio default risk. 
Taking the hazard rate to be small, the time horizon to be short and the defaults to be 
random, 

AT 

v^^\T)^{l-R) — , (43) 

and 



'''\T) /1-i? ^^^^ 



(i;(2)(T))3/2 y AT AT' 

where A is the actuarially fair spread. Hence with A held fixed the variance and skewness 
of the portfolio default loss probability distribution decrease as R increases. 



Since we are using the approximate formulas in equations (|4^) and (42) we would not 
get exact agreement with Table Ila even if i? = 0. For example with T = lyr, N = 10, 
7 = -4, ^ = lOObp, h = 0.02, i = 0.5 and i? = using equations (ji^) and (||) yields, 
a2 = 0.44, 03 = 0.33, 04 = 0.31 and 05 = 0.30. On the other hand the first row of Table 
Ila is, 02 = 0.42, 03 = 0.31, 04 = 0.29 and 05 = 0.29. 



6 Concluding Remarks 

Default correlations have an important impact on portfolio default risk. Given a default 
model the tail of the default loss probability distribution is difficult to compute, often 
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involving numerical simulation]^ of rare events. The first few moments of the default loss 
distribution are easier to calculate, and in the default model we used this computation 
involved some simple numerical integration. More significantly, the first few terms in the 
moment expansion have a familiar meaning. Investors are used to working with the classic 
mean, variance and skewness measures and have developed intuition for them and confidence 
in them. In this paper we studied the utility of the first few moments of the default loss 
probability distribution for the portfolio allocation problem. 

The default model we use assumes a multivariate normal Copula function. However, the 
corporate assets themselves are not necessarily normal or lognormal and their probability 
distribution can have fat tails. When there are correlations the portfolio default loss dis- 
tribution does not approach a trivial^ or normal probability distribution as the number of 
firms in the portfolio goes to infinity. Correlations dramatically decrease the effectiveness 
of increasing the number of firms, A^, in the portfolio to reduce portfolio default risk. 

We explored the convergence of the moment expansion for the optimal fraction of cor- 
porate assets Uopt- Our work indicates that, for Oopt less than unity, the convergence of the 
moment expansion is quite good. The convergence of the moment expansion gets poorer 
as aopt gets larger. We also examined how Uopt depends on the level of risk aversion, the 
recovery fraction and investment time horizon. 

The value of Oopt depends on the utility function used. In this paper we used utility 
functions with constant relative risk aversion, 1 — 7. It is possible to make other choices 
and while the quantitative results will be different for most practical purposes we expect 
that the general qualitative results should continue to hold. 

Our main conclusions from the examples in Sections 4 and 5 are: 

• In the presence of correlated defaults, the default loss distribution is not normal. 
Nevertheless we find that retaining just the first few moments of the default loss 
probability distribution yields a good estimate of the optimal fraction of corporates. 
As aopt gets smaller the convergence of the moment expansion improves. For the 
examples in sections 4 and 5 with Oopt near 40% retaining just the mean, variance and 
skewness of the default loss probability distribution determined Oopt with a precision 
better than about 15%. 

• For small numbers of corporate bonds increasing the number of firms decreases port- 
folio default risk. However, correlations cause this effect to saturate as the number 
of bonds increases and for the examples in sections 4 and 5 (with production factor 
or asset correlations 0.50 > > 0.25) the optimal fraction of corporates does not 
increase significantly if N is increased beyond 100 firms. 

• Increasing the recovery fractions decreases portfolio risk and this is true even when the 
probabilities of default increase so that the actuarially fair spreads remain constant. 

• With probabilities of default given by Pi(T) = 1 — exp{—hiT) and the hazard rates 
hi constant, continuously compounding of a constant excess risk premium does not 
exactly compensate for the increase in the default probability with time. The optimal 

^See, for example, DufBe and Singleton (1999) and Das, Fong and Geng (2001). 
^"^The trivial distribution has 51 = occuring with unit probability. 
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fraction of corporates decreased in going from a T = lyr to a T = 5yr investment 
horizon. 

The authors would like to thank their colleagues at PIMCO for enlightening discussions. 
MBW was supported in part by the Department of Energy under contract DE-FG03-92- 
ER40701. 
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